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ABSTRACT 


On  the  baali  of  a  statlatlcal  atudy  of  existing  experimental  data 
ralativa  to  ballistic  impact,  engineering  formulae  ware  developed  which 
will  predict  with  an  accuracy  of  90T, 

(a)  the  depth  of  penetration 

(b)  the  volume  of  the  crater 

in  semi- infinite  targets  of  common  metals  over  a  wide  range  of  impact 
velocities.  The  application  of  these  equations  in  predicting  the  size 
and  velocity  of  projectiles  impacting  orbital  structures  is  also  indicated. 

A  semi- rational  penetration  expression  was  developed  from  a  work- 
energy  consideration  which  suggests  that  the  nonrecoverable  target  com¬ 
pression  and  shear  strain  energies  may  account  for  most  of  the  kinetic 
energy  of  the  projectile.  This  dimensional  model  agrees  quite  well  with 
the  general  features  of,  and  is  subject  to,  the  same  limitations  as 
existing  theoretical  and  empirical  models. 
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-  Constanta  (Eq.  21) 

-  constants 

-  series  of  constants  as  02,  •  .  . 

"  dllatatlonai  wave  velocity 

-  projectile  diameter 

-  energy 

.  Impact  kinetic  energy  of  projectile 

-  tensile  elongation  (strain)  at  fracture  (%) 

-  target  resistance  force 

•  modulus  of  rigidity 

•  sheer  strain  at  fracture  (X) 

•  Brlnell  hardness 
«  bulk  modulus 

-  constant 

>  mass 

>  crater  depth  measured  from  original  target  surface 

-  pressure 

>  radius 

-  shear  strength 

-  target  temperature 

-  time 

-  ultimate  strength  (tensile) 

-  shock  velocity 


V 


volun* 


crater  volume  below  original  targat  aurfaea 

velocity,  Impact  velocity 

yield  point  or  yield  atrength 

shear  strain 

Poisson's  ratio 

mass  density 

shear  stress 

SUBSCRIPTS 


backsplash 

crater 

Interface  between  projectile  and  target 

projectile 

pressure 

denotes  material  strain,  as  £5  (strain  energy) 
target 

denotes  undisturbed  condition  (ahead  of  shock) 
denotes  compressed  condition  (behind  shock) 


used  with  constants 


INTRODUCTIOM 


With  the  advent  of  ballistic  missiles,  satellites,  and  space  vehi¬ 
cles  a  great  deal  of  attention  has  understandably  been  centered  about 
hypervelocity  terminal  ballistics.  The  velocity  of  meteorites  varies 
from  11.2  to  72  km. /sec.  At  a  velocity  of  50  km. /sec.,  for  example,  a 
gram  of  material  would  have  a  kinetic  energy  equal  to  the  potential 
energy  of  360  grams  of  TNT.  Therefore,  meteoric  material  presents  one 
of  the  greatest  natural  hazards  to  which  satellites  and  space  vehicles 
may  be  subjected.  In  addition,  artificial  "meteorites"  with  high  rela¬ 
tive  kinetic  energies  appear  attractive  as  anti  ICBM  defensive  weapons. 

The  highest  velocity  discreet  particles  so  far  reported  in  terminal 
ballistic  studies  are  9  km. /sec.  fragments  from  shaped  charges.  Not  only 
does  a  large  gap  exist  between  present  experimental  techniques  and  the 
maximum  velocity  of  meteoric  material,  but  the  terminal  ballistics  of 
projectiles  traveling  less  than  9  km. /sec.  are  poorly  understood. 

The  purpose  of  this  study  is  to  gather  and  assemble  existing  data  on 
ballistic  Impact  and  on  material  failure,  especially  at  high  impact  veloc¬ 
ities  or  large  loading  -  to  establish  the  relative  importance  of  such 
factors  as  projectile  velocity,  mass,  projectile-target  strengths,  duc¬ 
tilities,  densities,  compressibilities,  etc.,  and  to  use  this  information 
to  deduce  the  approximate  mathematical  relationships  of  critical  factors 
as  the  target  structure  responds  to  impact  and  is  penetrated. 
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GENERAL  ASPECTS  OF  TARGET  BEHAVIOR 


From  an  engineering  point  of  view,  the  behavior  of  a  target  under 
balllatlc  Impact  might  be  aeparated  Into  three  categorlea  ...  penetration, 
penetratlon-plua-scabblng,  and  perforation.  Scabbing,  of  courae,  la  luat 
one  phenomenon  of  tenalle  fracture  reaultlng  from  atresa  wave  Interference, 

In  addition  to  tenalon  (or  compreaalon)  the  targeta  are  aubject  to  ahear, 

but  In  each  caae  the  reaultlng  fracture  can  uaually  be  attributed  to  tran- 

alent  streaa  condltlona  that  develop  aa  a  reault  of  atresa  wave  Interferencea.  * 

Some  materlala  (ateela)  exhibit  a  time  delay  before  the  Initiation  of  plaa- 

tic  flow  which  tenda  to  produce  a  more  brittle  type  of  fracture,  while 

other  metala  (aluminum,  braaa)  do  not  ahow  theae  time- dependent  plaatic  p 

flow  propertlea.  Although  there  are  a  multiplicity  of  factora  which  may 

affect  target  behavior,  the  factora  that  determine  which  of  the  three 

behavior  categorlea  will  reault  for  a  given  aet  of  Impact  condltlona  and 

materials  are  principally  the  target  thickness  and  whether  the  target 

material  behaves  substantially  aa  a  brittle  or  ductile  material. 

The  problem  of  predicting  material  resistance  to  high  speed  Impact 
la  difficult  largely  becauae  the  question  of  ductility  versus  brittleness 
la  a  relative  and  variable  one.  While  some  materials  are  more  brittle  or 
ductile  than  others  for  a  given  set  of  conditions,  a  particular  material 
may  behave  In  a  brittle  or  ductile  fashion  depending  on  such  conditions 
as  temperature,  strain  rate,  and  pressure.  In  general,  a  particular 
material  tends  to  become  less  ductile  at  higher  strain  rates  and  lower 
temperatures  and  more  ductile  under  high  pressures.  There  la,  at  least 
for  ferrous  metals,  a  definite  transition  temperature  range  separating 
brittle  and  ductile  behavior.  Rapid  strain  rates  tend  to  raise  this 
transition  temperature  range  so  that  decidedly  brittle  behavior  occurs 
at  higher  temperatures  when  subjected  to  higher  strain  rates.  This  transi¬ 
tion  range  la  not  generally  found  among  non-ferrous  metala. 


DISCUSSION  OF  CRATER  PARAMETERS 


Metal  targets  tend  to  pass  from  "penetration"  to  "penetratlon-plua- 
scabblng"  to  "perforation"  types  of  behavior  as  the  target  thickness 
decreases.  However,  the  mechanism  of  these  different  behaviors  and  their 
transition  from  one  to  the  other  are  quite  different  in  brittle  and  ductile 
targets.  A  relatively  thick  ductile  target  will.  In  general,  undergo  a 
more  or  less  hemispherical  crater  formation  under  hypervelocity  ballistic 
impact.  In  thla  case,  the  depth  of  penetration  and  the  crater  volume  would 
assume  equal  status  as  the  design  variable  dependent  on  several  different 
Impact  and  material  parameters.  The  mechanism  of  crater  formation  In  a 
thick  brittle  target  Is  somewhat  different  in  that  the  brittle  target  may 
tend  to  spall  In  the  vicinity  of  the  crater  rim  and  Is  more  susceptible 
to  scabbing  fracture.  Since  the  energy  necessary  to  propagate  a  crack 
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It  quite  snail,  It  teems  reasonable  that  spall  may  have  a  primary  effect 
on  crater  volume,  but  only  a  secondary  effect  on  penetration.  An  example 
of  this  secondary  effect  on  penetration  would  be  that  the  spalled  away 
material  might  give  a  larger  solid  angle  for  the  low  resistance  ejection 
of  material  from  the  growing  crater.  As  long  as  the  Impact  velocity  Is 
high  enough.  It  might  appear  that  depth  of  penetration  Into  thick  tar¬ 
gets  would  be  less  sensitive  to  the  ductility  of  the  target  than  would 
crater  volume. 


DISCUSSION  OF  MATERIAL  PARAMETERS 

For  given  projectile  and  target  geometry  and  Impact  conditions, 
target  behavior  la  a  function  of  the  pertinent  material  properties  that 
are  operative  under  the  prevailing  Impact  conditions,  which  Include  high 
pressures,  high  strain  rates  and  possible  high  or  low  ambient  target 
temperatures.  Since  the  fundamental  material 'properties  used  In  this 
study,  of  necessity,  refer  to  conditions  which  are  greatly  different 
from  those  under  hypervelocity  Impact  loading,  some  justification,  or 
at  least  discussion,  of  the  adequacy  of  their  use  Is  required. 

The  projectile  and  target  material  properties  used  In  the  study 
Include  densities  under  normal  pressures  and  temperatures,  and  mechanical 
properties  under  static  or  quasl-statlc  loading.  With  the  exception  of 
shear  strength  and  Brlnell  hardness,  all  of  the  mechanical  properties  used 
are  tensile  properties. 

While  It  Is  well  known  that  material  behavior  changes  radically 
under  conditions  of  high  pressure,  high  strain  rate  and  large  tempera¬ 
ture  changes,  it  Is  believed  that  the  properties  used  In  this  study  are, 
generally  speaking.  Indicative  of  the  behavior  under  conditions  of  high 
speed  Impact.  That  Is,  higher  strain  rates  generally  produce  greater 
material  strengths,  less  compressibility  and  less  ductility.  Tempera¬ 
ture  decreases.  In  general,  have  similar  effects  on  material  strength, 
compressibility  and  ductility.  Also  decreasing  temperatures  generally 
Increase  material  hardness.  With  regard  to  high  and  low  target  tempera¬ 
tures,  It  appears  that  high  ambient  target  temperatures  (below  the 
melting  point  of  the  material)  are  probably  not  as  Important  as  low 
temperature  effects,  since  material  behavior  at  higher  temperatures 
usually  takes  the  form  of  a  delayed  creep  or  flow,  which  would  not  seem 
to  be  pronounced  under  very  rapid  transient  loads.  High  pressures 
generally  Increase  material  strength,  ductility  and  hardness,  and  de¬ 
crease  compressibility.  The  use.  In  a  comparative  manner,  of  the  more 
accurately  defined  static  and  low  strain  rate  material  properties  seems 
to  be  further  justified  In  light  of  the  knowledge  that  many  factors  such 
as  thermal  effects  and  modification  of  the  crystalline  structure  of  the 
target  material  by  transient  stress  waves  may  have  a  decided  influence 
on  behavior,  apart  from  the  mechanical  properties  of  the  target  material, 
and  factors  such  as  these  are  little  understood  or  difficult  to  evaluate. 


Th«  following  it  •  tuiiinary  of  tho  gonoral  offoets  of  InerMtod  otraln 
ratos  and  praaaurai  and  daeraatad  tanparaturos  on  oatarlal  bahavlor: 


Strength  (Generally) 
Conipreeslblllty  (Generally) 
Ductility  (Generally) 
Brlnell  Hardneaa 


Increased  Strain  Kata 


Greater 

Less 

Lwaa 


Decreaeed 

Temp. 

Greater 

Less 

Less 

Greater 


Increased 

Press. 

Greater 

Less 

Greater 

Greater 


The  Brlnell  definition  of  hardness  Is  specified  since  hardness  Is  defined 
in  various  ways,  some  of  which  may  not  follow  this  general  behavior  (for 
example,  when  defined  ae  energy  absorbing  capacity). 

The  previous  discussion  has  sunanarleed  a  few  of  the  more  pertinent 
aspects  of  material  behavior  under  conditions  which,  for  low  tempera¬ 
tures,  probably  approximate  those  found  In  the  hypervelocity  Impact  of 
a  meteoric  particle  with  a  apace  vehicle.  However,  our  experimental 
Information  on  the  effect  of  high  strain  rates  and  high  pressures  on 
fundamental  material  behavior  falls  far  short  of  those  experienced  In 
hypervelocity  Impact  loading.  Also,  It  would  be  foolhardy  not  to  expect 
that  some  contradictions  to  the  above  general  rules  exist.  For  example, 
certain  grades  of  mild  steel  are  known  to  have  a  so-called  "brittle  range" 
under  static  loading  (ordinary  room  temperature  to  about  500  degrees  Fah¬ 
renheit)  where  an  Increase  In  temperature  Is  accompanied  by  an  Increase 
In  strength  and  brittleness.  This,  however.  Is  abnormsl  when  compared 
with  other  temperature  ranges  for  the  same  material  and  Is  not  found  at 
all  In  most  other  metals. 

• 

Bjork^  has  pointed  out  an  apparent  Inconsistency  In  relative  material 
compressibility  under  high  pressures  as  compared  with  low  pressures,  al¬ 
though  It  does  not  violate  the  general  statement  that  compressibility 
decreases  with  pressure.  He  showed  that  the  hugonlots  for  some  metals 
tend  to  "cross  over".  Specifically,  the  graphs  of  pressure  versus  rela¬ 
tive  density  (/^^,  )  for  lead  and  aluminum  cross.  At  low  pressures 
aluminum  Is  more  resistant  to  density  change  than  lead,  but  at  high 
pressures  this  Is  reversed.  While  a  theoretical  study  of  the  mechanics 
of  Impact  behavior  must  necessarily  consider  the  hugonlots.  It  Is  probable 
that  in  a  statistical  study  this  effect  can  be  assumed  to  be  reflected 
generally  In  the  relative  material  densities,  strengths,  compressibilities, 
etc.  The  Important  fact  In  this  regard  seems  to  be  that  the  hugonlots  for 
different  materials  have  the  same  general  shape  and  that  eonpresilbillty 
generally  decreases  with  Increased  preasure. 


agntlCAL  MODEL 


Approxlmataly  1800  exp«rlm«ntal  data  polnta  have  bean  accwmilated 
from  the  original  papers  of  many  Imres tiga tors  and  analyzed  during  the 
course  of  this  atddy.  Moat  of  the  accumulated  data  (with  the  major 
exception  of  the  higher  velocity  shots,  especially  In  aluminum)  have 
been  compiled  and  analysed  by  BruceSl  or  Herrmann  and  Jonea^.  Of  the 
1800  shots  only  about  half  have  proven  really  useful  and  there  are 
many  discrepancies  among  these  useful  shots.  The  majority  of  shots 
were  eliminated  for  two  reasons:  (1)  failure  of  the  Investigator  to 
adequately  specify  the  material  of  the  target,  and  (2)  the  large 
divergence  of  some  Impact  data  from  other  data  taken  by  other  Investi¬ 
gators  Into  the  same  materials  under  very  similar  Impact  conditions. 
Many  experlmentors  have  chosen  to  present  their  data  only  on  curves 
and  graphs  rather  than  In  tables  of  raw  data,  thus  Increasing  the 
chance  of  recording  and  reading  errors. 

It  was  decided  to  relate  the  dependent  variables  (penetration 
depth  and  crater  volume  in  semi- infinite*  targets)  to  the  Independent 
variables  through  a  simple  power  law.  Ho  additional  assumptions  were 
made  regarding  the  process  of  ballistic  Impact  or  the  shape  of  the 
craters  formed.  The  analysis  has  produced  engineering  equations  which 
will  predict,  within  the  limits  of  the  data,  the  penetration  depth  and 
crater  volume  in  common  metallic  targets  with  an  accuracy  of  about  901. 
By  the  very  manner  In  which  experimental  data  must  be  gathered  and 
crater  parameters  must  be  measured.  It  Is  doubtful  that  the  reproduci¬ 
bility  of  the  original  data  Is  greater  than  901. 


FRELIMIH^RY  EHPIRICAL  ANALYSIS  OF  Pc 


A  preliminary  analysis  vis  performed  on  data  compiled  prior  to 
December,  1961  was  reported  In  the  Second^  and  Thlrd^  Qviarterly  Reports. 
Since  the  distribution  of  target  materials  as  a  function  of  velocity 
was  quite  biased,  the  detailed  results  of  this  preliminary  analysis 
are  not  very  useful.  It  was  performed  more  to  establish  methods  for 
later  analyses  and  the  detailed  results  will  not  be  reproduced  In  this 
raport. 


Penetration  depth  was  considered  a  function  of  ten  Independent 
variables  according  to  the  equation 


Pc 


ko  V 


ki  kj  k3  k4 
t  /Op  Vp 


l^eml-lnflnlte  according  to  the  rule-of-thumb  published  by  Klnard, 
at  al',  that  the  penetration  Is  no  greater  than  20%  of  the  target  thickness 


5 


This  •quatlon  was  reduced  to  Its  eeeoclated  linear  form 

log  Pc  ■  log  ko  +  ki  log  V  + . +  kjQ  log 

and  the  method  of  least  squares  was  used  to  determlns  a  "best"  set  of 
values  for  the  variable  exponents  (coefficients)  ki  based  on  minimising 
ths  value  of 


2  log  (Pc)experimental 
1*0 


log  (Pc) calculated 
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The  Independent  variables  were  dropped  one-by>one  starting  with  Dp  and 
the  effect  on  the  kj^  was  noted. 

After  the  Icj  were  computed,  various  statistical  teats  were  made  to 
determine  the  validity  of  the  results  obtained.  The  per  cent  of  the 
variance  in  the  log  of  dependent  variable  (Pc)  that  is  explained  by  the 
logs  of  the  Independent  variables  is  given  by  the  multiple  correlation 
coefficient.  These  multiple  correlation  coefficients  were  calculated 
after  each  independent  variable  was  dropped.  An  analysis  of  this  vari> 
ance, using  the  F  distribution^ tests  the  null  hypothesis  that  oil  /3  i  are 
sero.  Separate  tests  of  the  individual  hypothesis  /3.  ■  0  were  made 
using  the  t  distribution.  Variables  corresponding  to^the/di  that  were 
not  significant  by  the  preceding  test  were  then  removed  and  separate  tests 
were  made  on  each  of  the  removed  variables  to  see  whether  the  adding  of 
this  variable  improved  the  fit.  Individual  correlation  coefficients 
between  all  Independent  variables  were  also  calculated. 


The  significance  of  the  target  yield  strength  above  many  of  the  other 
variables  contributing  to  the  multiple  correlation  coefficients  was  noted. 
However,  the  target  yield  strength  (Yt),  the  dllatational  wave  velocity  in 
the  target  (c^),  and  the  target  ultimate  strength  (U^)  are  statistically 
indistinguishable  as  evidenced  by  the  fact  that  the  calculated  individual 
correlation  between  them  is  over  95X.  The  diameter  of  the  projectile  (Dp) 
is  also  statistically  indistinguishable  from  the  volume  of  the  projectile 
(Vp)  so  that  the  effect  of  dropping  (Dp)  was  insignificant.  With  the 
exception  of  the  above  correlations  ana  the  dependence  of  ct  on  /9t> 
rest  of  the  individual  correlation  coefficients  showed  a  low  degree  of 
interdependence . 


The  preliminary  analysis  was  made  on  approximately  1200  experimental 
shots  only  half  of  which  were  reported  in  enough  detail  to  enable  assign¬ 
ment  of  strength  parameters  to  the  target  and  projectile.  Of  the  resulting 
600  shots,  only  174  were  performed  at  velocities  greater  than  the  bulk  wave 
velocity  in  the  target  medium  and  over  90%  of  these  "supersonic"  shots  were 
in  lead  targets  where  the  bulk  wave  velocity  is  quite  low  (about  5400  feet 
per  second).  This  analysis,  then,  seemed  to  indicate  that  ths  crater  depth 
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could  be  adequately  predicted  on  the  baeis  of  4  or  5  out  of  the  ten  inde« 
pendent  varlablea  choaan,  at  least  for  Impact  valocltles  below  the  bulk 
or  dtlatatlonal  wave  velocity  In  the  target  wedlum. 


ANALYSIS  OF  Pe/Dn 

An  analysis  of  the  penetration  depth  In  seml-lnflnlte  targets  divided 
by  the  diameter  of  the  projectile  (Pc/Dp)  vas  performed  on  1,270  shots 
which  were  reported  in  enough  detail  to  assign  strength  parameters  to  the 
target  materials.  The  experimental  data  was  divided  Into  two  parts: 

(1)  those  having  impact  velocities  below  V3  (Hopkins  and  Kolsl^),  end 

(2)  those  with  impeet  velocities  near  or  eMve  V3,  where 

*3  -  y 

The  distribution  of  target  materials  as  a  function  of  Impact  velocity  Is 
given  In  Figure  1*  for  the  "low"  velocity  set  (985  shots) .  Since  a  larga 
number  of  recording  errors  were  discovered  In  the  high  velocity  data,  tha 
reaults  of  this  run  (recorded  In  the  Fourth  Quarterly  Report^)  have  been 
eliminated  from  the  present  report. 

As  was  done  In  the  previous  analysis  of  (Pc)»  (Pc/^)  fitted  to 
a  sisple  power  law  formula 


k^  k2  kj 

^c/Dp-ko^  /^t  . 

with  no  additional  assumptions  being  made.  This  aquation  was  then  reduced 
to  Its  associated  linear  form 


^c/Dp  ■  ^'o  *^1  V  +  k2  log/Oj.  +, 


and  the  method  of  least  squares  was  used  to  determine  a  "best"  set  of  values 
for  the  variable  exponents  (coefficients)  ki  based  on  minimising  the  value 


2  F* 

t«0  __ 


c /Dp^ exper  Ime ntal 


The  Independent  variables  were  dropped  one-by-one,  starting  with  the  last, 
and  the  effect  on  the  resmlnlng  ki  determined.  These  are  summarised  for 


*  Appendix  B 
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the  low  velocity  dete  in  Tables  1,  II,  III  end  IV,  which  eppeer  la 
Appendix  B. 

After  the  ki  were  computed,  the  usual  statistical  tests  were  made 
to  determine  the  validity  of  the  results.  The  per  cent  of  variance  In 
log  Pc/Dp  that  is  explained  by  the  variation  In  the  logs  of  the  indepen¬ 
dent  variables  is  given  by  the  multiple  correlation  coefficients  listed 
In  each  table.  Individual  correlation  coefficients  between  all  variables 
were  also  computed  on  the  low  velocity  shots.  These  results  appear  as 
Table  V,  also  In  Appendix  B. 

From  Table  V,  it  may  be  seen  that  the  close  correlation  between  the 
dilatational  wave  velocity  (ct),  the  ultimate  strength  (Ut),  the  shear 
strength  (St),  the  Brinell  hardness  (Ht),  and  the  yield  strength  (T^) 
render  these  independent  variables  statistically  indistinguishable.  In 
addition,  the  correlation  coefficient  between  the  bulk  siodulus  (Kt)  and 
each  of  the  above  variables  is  greater  than  0.50.  Any  of  the  empirical 
power  law  equations  containing  more  than  one  of  theae  closely  correlated 
variables  is  likely  to  be  completely  misleading  in  the  coefficients 
(exponents)  of  the  correlated  variables.  All  the  rest  of  the  variables 
Including  the  per  cent  elongation  (e^)  were  statistically  Independant. 


ANALYSIS  OF  ?c/V^ 

For  this  analysis,  the  cube  root  of  the  projectile  volume  was  used 
to  render  the  penetration  non-dimensional.  This  not  only  makes  the 
results  of  the  empirical  analysis  assume  the  same  form  as  the  results 
of  the  Dimensional  Model  which  follows,  but  also  eliminates  the  ambi¬ 
guity  involved  in  defining  the  projectile  diameter  when  the  projectile 
is  not  exactly  spherical. 

Before  this  analysis  was  run,  elaborate  steps  were  taken  to  "clean 
up"  the  accumulated  experimental  data.  First,  the  "high"  velocity  data 
were  completely  rechecked  for  recording  errors  and  shots  at  other  than 
normal  impact  angles  were  eliminated.  The  body  of  "low"  velocity  data, 
being  much  larger,  was  checked  on  the  computer.  The  equation  from  Table  III, 

-3  1.06  -.365  .984  .400  .586  .010  .035  .252 

Pc,  -  8.43  X  10  V  /o  t  ^  St  Tt  et  Vp  Ct 

Up 

which  yielded  a  multiple  correlation  coefficient  of  0.846  on  the  previous 
analysis,  was  used  to  calculate  the  depth  of  penetration  from  the  independ¬ 
ent  variables  of  each  of  the  "low"  velocity  shots.  The  calculated  penetra¬ 
tion  depth  was  then  cosipared  with  the  experimental  penetration  depth  and  a 
per  cant  error  was  computed  for  each  shot.  All  data  points  with  an  error 
of  greater  than  50%  were  then  scrutinized  for  recording  or  computer  errors. 
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Shots  which  showed  a  large  deviation  and  whose  values  had  been  read  from 
curves  rather  than  tables  In  the  original  papers  were  Inmed lately  eliminated. 
Other  shots,  which  did  not  seem  to  fit  with  the  rest  of  the  experimental 
data  taken  from  their  particular  source  were  also  eliminated.  Many  shots 
with  errors  greater  than  50%  were  deemed  correct  and  were  left  In.  Finally, 
the  small  number  of  shots  In  which  the  target  temperature  or  the  Impact 
angle  were  other  than  ambient  and  90°,  respectively,  were  eliminated. 

After  "cleaning",  a  total  of  1008  shots  remained  .  742  classified 

as  "low"  velocity  and  266  classified  as  "high"  velocity.  The  distributions 
of  target  materials  as  a  function  of  Impact  velocity  for  each  of  these 
compilations  are  shown  In  Figures  2  and  3  of  Appendix  B. 

The  "low"  velocity  and  the  "high"  velocity  data  wera  Individually  fit 
to  the  simple  power  law  formula  In  the  form 

ki  kj  k3 

Pc/'^  -kov  /Ot  /Op  . 

P 

with  the  digital  computer  programmed  to  calculate  least  squares  values 
of  the  ki  and  the  multiple  correlation  coefficients  between  the  dependent 
variable  and  all  of  the  Independent  variables  taken  together.  As  was  done 
In  the  previous  runs,  the  kj^  were  dropped  one  by  one  and  the  effect  on  the 
multiple  correlation  coefficient  was  noted.  The  product  of  the  shear 
strength  and  the  per  cent  elongation  (Se)  was  used  as  one  of  the  variables 
In  this  analysis  because  of  Its  appearance  In  the  Dimensional  Model  which 

follows.  The  results  of  eight  separate  runs  appear  In  Appendix  B  . . 

Tables  VI,  VII,  VIII,  and  IX  for  the  "low"  velocity  data  and  Tables  X,  XI, 
XII,  and  XIII  for  the  "high"  velocity  data. 

The  multiple  correlation  coefficients  In  Tables  1-8  show  that  using 
four  Independent  variables  (l.e.,  v, /9p,  S^,  and /o  ()  It  Is  possible  to 
explain  92%  of  the  variation  In  penetration  depth  using  the  low  velocity 
data  and  87%  of  the  variation  In  penetration  depth  using  the  high  velocity 
data. 

ANALYSIS  OF  Vc/Ep 

Finally,  the  ratio  of  crater  volume  to  kinetic  energy  was  used  as  the 
dependent  variable  In  the  equation 

k.  k  k,  k 
Vc/e  -  ^  /®p  St  /Ot 

P 

As  usual,  the  were  found  by  the  method  of  least  squares.  Since  It  had 
been  well  established  by  this  time  that  only  four  Independent  variables 
(Including  a  strength  parameter)  are  necessary  to  adequately  describe  the 
penetration,  the  variables  were  not  dropped  one  by  one.  The  "low"  (Fig.  2) 
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and  "high"  (Fig.  3)  velocity  data  were  substituted  separately.  Because 
the  crater  volume  was  not  measured  on  all  data  points  for  which  the  pene¬ 
tration  depth  was  obtained,  only  543  "low"  velocity  shots  and  220  "high" 
velocity  shots  contributed  to  this  analysis.  The  results  appear  In 
Table  XV  of  Appendix  B. 

QUALITATIVE  DISCUSSION  OF  THE  EMPIRICAL  MODEL 

Reference  to  Figures  2  and  3  In  Appendix  B  shows  that  the  "low" 
velocity  shots  are  fairly  well  distributed  among  target  materials.  Each 
bar  on  Figure  2  represents  the  total  number  of  shots  over  a  velocity  range 
of  one  thousand  feet  per  second  while  each  bar  on  Figure  3  represents  a 
range  of  two  thousand  feet  per  second.  The  relatively  large  manber  of 
shots  Into  lead,  copper,  and  2024  T4  aluminum  represent  soft,  medium, 
and  fairly  hard  materials,  respectively.  Figure  3  shows  a  very  large 
peak  from  1000  to  5000  feet  per  second,  but  any  attempt  to  flatten  the 
distribution  as  a  function  of  velocity  would  have  meant  the  elimination 
of  a  tremendous  number  of  shots  with  a  consequent  decrease  In  the  statis¬ 
tical  confidence.  Figure  3  points  out  a  definite  need  for  experimental 
data  over  the  velocity  range  of  17,000  to  27,000  feet  per  second.  Light 
gas  guns  have  been  developed  with  muzzle  velocities  up  to  25,000  feet 
per  second,  but  so  far  no  cratering  data  has  appeared  In  the  literature 
at  these  high  velocities.  The  experiments  conducted  at  velocities  above 
25,  000  feet  per  second  all  employed  shaped  charge  techniques  (Aerojet 
General  for  the  Terminal  Ballistic  Division  at  Eglin  AFB,  Florida). It 
would  be  nice  to  see  more  impacts  into  other  than  2023  T4  targets. 
Finally,  there  is  a  need  for  data  at  much  higher  velocities. 

In  practically  all  of  the  data  used  In  this  empirical  study,  the 
projectile  shape  Is  spherical  or  near-spherical  (l.e.^  cylinder  with 
1/D  -  1  or  fragment) .  At  lower  Impact  velocities  thejre  Is  no  doubt  that 
the  shape  of  the  projectile  has  a  pronounced  Influence  on  cratering. 

As  the  Impact  velocity  Increases,  It  would  seem  reasonable  to  assume 
that  the  shape  of  the  projectile  becomes  less  and  less  Important.  How¬ 
ever,  over  the  velocity  range  considered  In  this  study  one  would  still 
expect  the  length  to  diameter  ratio  (t/D)  of  the  projectile  to  have  some 
effect.  Until  some  high  Impact  velocity  (possibly  equal  to  the  sum  of 
the  dllatstlonal  wave  velocity  In  the  projectile  and  target  materials) 

Is  reached.  It  will  be  possible  for  the  tall  of  a  projectile  to  be  aware 
that  something  has  happened  to  the  head.  Little  experimental  work  has 
been  done  with  high  &/D  ratios  and  It  Is  not  likely  that  any  will  appear 
for  very  high  velocities.  It  is  difficult  enough  to  obtain  a  discreet 
projectile  at  very  high  velocities,  let  alone  control  Its  shape.  Shaped 
charge  jets  might  approximate  high  |/D  ratio  projectiles  were  It  not  for 
the  fact  that  the  jet  Is  moving  with  nonuniform  velocity  along  Its  length. 
Fortunately,  from  a  practical  standpoint.  It  Is  quite  unlikely  that  space 
debris  will  appear  In  the  form  of  long  needles. 

The  fact  that  the  multiple  correlation  coefficients  for  the  "high" 
velocity  data  (Tables  X  through  XIII)  are  somewhat  lower  than  those  for 
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th*  "low"  velocity  date  (Tables  VI  through  IX)  might  have  been  anticipated. 
Higher  impact  velocities  become  increasingly  more  difficult  to  measure 
accurately.  In  addition,  the  mass  of  the  projectile  in  the  shaped  charge 
experiments  (all  above  25,000  feet  per  second)  was  determined  from  orthog¬ 
onal  photographs.  For  these  and  other  reasons  the  higher  velocity  data 
is  less  reproducible  than  the  lower  velocity  data. 

It  might  be  inferred  from  the  results  of  this  study  that  as  long  as 
four  independent  variables  (  v  •  /Op>  /Ot>  plus  one  "strength"  parameter) 
are  used  in  the  prediction  formula,  increasing  the  number  of  independent 
variables  has  little  effect  on  the  correlation.  This  is  not  true.  As 
pointed  out  earlier  in  connection  with  Table  V,  all  of  the  mechanical 
properties  used  in  this  study  seem  to  have  high  correlation  with  one 
another.  One  would  not  expect  the  addition  of  variables  highly  corre¬ 
lated  with  existing  variables  to  Improve  the  multiple  correlation.  Tar¬ 
get  temperature  was  used  as  a  variable  in  the  analysis  of  and  P^/D., 
but  the  number  of  shots  in  which  the  target  temperature  varied  was  so*^ 
small,  relative  to  the  total  number  of  shots,  that  dropping  the  tempera¬ 
ture  as  a  variable  had  little  effect  on  the  multiple  correlation  coeffi¬ 
cient.  Likewise,  impact  angle  was  not  used  as  a  variable  due  to  the  small 
number  of  shots  at  impact  angles  other  than  90°. 

It  is  Interesting  to  note  that  the  "low"  velocity  data  indicate  a 
dependence  of  cratering  on  the  first  power  of  the  impact  velocity.  This 
linear  dependence  of  (P^)  on  (v)  agrees  with  the  diacusaion  of  the  simple 
power  law 

n 

P^  ■  k  V 
c 

Q 

by  Herrmann  and  Jones  in  which  they  find  that  the  experimental  data,  at 
least  for  aluminum,  fits  n  »  1  over  the  approximate  velocity  of  7500  to 
12,000  feet  per  second.  They  also  point  out  that  in  the  velocity  range 
of  10,000  to  20,000  feet  per  second  penetration  may  be  approximated  by 
the  above  equation  with  n  -  2/3,  and  from  30,000  to  200,000  it  may  be 
approximated  with  n  -  1/3.  This  type  of  high  velocity  behavior  seems  to 
be  supported  by  the  fact  that  the  "high"  velocity  data  (Tables  X  through 
XIII)  indicate  a  dependence  of  (P^/Vp*^!  )  on  a  power  of  (v)  somewhat  less 
than  1/2.  It  might  also  be  said  that  the  decreasing  power  of  /Q  as 
velocity  increases  indicates  that  the  data  is  approaching,  but  hSs  not 
reached,  the  velocity  regime  in  which  Bjork's  hydrodynamic  theory  subtends. 

The  low  velocity  data  also  seem  to  indicate  a  dependence  of  (Pj./Dp) 
or  (Pg/Vp'^  )  on  the  first  power  of  the  projectile  density  (/9p).  This, 
coupled  with  the  linear  dependence  on  velocity  mentioned  abwe,'^fits 
Herrmann  and  Jones'  assumption  of  a  resiatlve  force  dependent  on  velocity  as 

P  -  k  Dp^  V 
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which  iMdc  to  «  penetration  expraaelon  of  the  fora 


where  k  Is  a  constant.  As  they  point  out,  an  equation  somewhat  like 
the  above  can  be  used  to  fit  low  velocity  Impact  data. 

The  high  velocity  data  show  a  dependence  of  )  on  to  a 

power  of  about  2/3.  This  la  a  definite  disagreement  with  the  widely  used 
empirical  expression 


P 


■  k  v*^* 


which  Is  obtained  from  the  observation^^*  ^^*  ^^*  ^  ^®*  ^*  ® 

that  the  crater  volume  Is  proportional  to  the  kinetic  energy,  coupled  with 
the  observation  ^^*  that  the  craters  formed  are  hemispher¬ 

ical.  Palmer  points  out  that  all  Utah  cratering  studies  seem  to  best  fit 
an  equation  similar  to  the  above  equation  In  which  the  exponent  of  Is  1/2. 

Both  the  "high”  and  "low"  velocity  data  point  to  a  real  dependence 
of  cratering  on  at  least  one  target  "strength"  parameter.  The  "low" 
velocity  data  (Tables  VI,  VIII,  and  IX)  favor  an  Inverse  dependence  of 
(Pc/Vp'^i  )  on  S(.,  Hf,  or  (Se)j^  to  a  power  slightly  more  than  1/3  while 
the  "nigh"  velocity  data  (Tables  X,  XU,  XIII)  favor  an  Inverse  power 
slightly  less  than  1/3.  This  Is  In  excellent  agreement  with  Palmer's^^ 
observation  that  the  crater  depth  Is  Inversely  proportional  to  1/3  power 
of  the  target  shear  strength,  and  Herrmann  and  Jones^  and  Elchelberger's^^ 
observations  that  It  Is  Inversely  proportional  to  the  1/3  power  of  the 
target  Brlnell  hardness.  Table  IV  favors  an  Inverse  dependence  on  Kt 
to  a  power  slightly  greater  than  1/2  when  It  Is  used  as  the  target  material 
parameter. 

UANTITATIVE  USE  OF  THE  EMPIRICAL  MODEL 


The  primary  aim  of  this  study  was  to  provide  at  least  first  approxi¬ 
mation  formulae  by  which  penetration  Into  various  materials  might  be  pre¬ 
dicted.  Certainly  an  accuracy  of  90%  Is  better  than  a  first  approximation. 

The  tables  of  data  are  all  collected  In  Appendix  B  to  facilitate  their 
use.  Each  group  of  tables  Is  preceded  by  the  appropriate  distribution  of 
target  materials  as  a  function  of  velocity  to  Indicate  generally  the  limits 
of  applicability.  It  should  be  pointed  out  that,  although  the  dependent 
variable  Is  expressed  In  nondlmenslonal  form,  the  Independent  variables 
are  not  In  nondlmenslonal  groupings.  Therefore,  the  constants  k^  will  have 
odd  sets  of  units  which  depend  upon  the  Independent  variables  used  and  thalr 
specific  units.  For  this  reason.  Table  XIV  Is  Included  to  give  the  material 
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parameters  used  In  this  study  in  their  proper  units  ilor  insertion  into 
any  of  the  empirical  formulae. 

For  example,  the  "low"  velocity  relationship  yielding  the  highest 
correlation  becomes  (Table  VI,  line  4) 


when  only  four  independent  variables  are  used.  In  order  to  obtain  a 
dimensionless  number  for  the  ratio  of  P^/V  ,  it  is  necessary  to  sub¬ 
stitute  the  impact  velocity  in  feet  per  second,  the  densities  in  pounds 
per  cubic  inch,  and  the  shear  strength  in  pounds  per  square  inch. 


Taking  one  shot  from  the  Utah  data  (Turner,  Palmer,  and  Grow,  Uni¬ 
versity  of  Utah  Technical  Report  UU-5.  August  2,  1960,  page  77,  line  6) 
where  the  target  material  is  lead  and  the  projectile  material  is  aluminum. 
The  experimental  data  gives 

Pc  -  0.154  in.  Vc  -  .011  in.^ 

V  -  3150  ft. /sec.  Vp  -  .00347  In.^ 

From  Table  XIV 

/O  -  .1001  Ib./ln.^  S..  -  1800  lb./ln.2 
'  P  t 


/O^  -  .4097  lb. /in. 3 

Substituting  these  independent  variables  into  the  equation  yields 

.893  .979  -.457  -.350 

Pg/(. 00347)  -  .172  (3150)  (.1001)  (1800)  (.4097) 


so  that  Pg  ■  0.161  in. 

which  represents  a  4.57L  error  when  compared  with  the  experimental  penetra¬ 
tion  depth. 

Using  the  second  line  of  Table  XV 

-.143  .374  -.814  -.227 

0216  V 

2 

and  the  fact  that  Ep  ■■  1/2  Vp  v  ,  the  equation  for  crater  volume  be¬ 
comes 
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1.857  1.374  -.814  -.227 

V^,  -  .0108  V 

''P 

Substitution  of  ths  independent  variables  above  yields 

Vc  -  0.013  ln.3 

which  represents  a  20%  error  when  compared  with  the  experimental  crater 
volume . 

It  is  Just  coincidental  that  a  shot  was  chosen  whose  crater  depth 
was  predicted  within  the  multiple  correlation  coefficient  and  whose  i 

crater  volume  was  not. 

It  is  interesting  to  note  that  the  equations  developed  in  this  report 
may  be  used  to  analyze  cratering  in  materials  recovered  from  orbit,  and 
thereby  make  predictions  on  the  velocity  and  size  of  the  impacting  projec¬ 
tile.  If  the  projectile  density  and  the  target  material  are  known,  and 
if  the  penetration  depth  and  crater  volume  are  measured  in  the  recovered 
structure,  it  is  possible  to  predict  the  impact  velocity  and  projectile 
volume.  This  can  be  done  with  any  degree  of  confidence  only  within  the 
impact  velocity  and  material  density  ranges  covered  by  the  experimental 
data. 


The  empirical  equations  succeed  in  cutting  across  Impact  and  material 
parameters  only  as  well  as  the  experimental  data  has  been  distributed 
across  these  same  parameters.  However,  the  high  degree  of  predictability 
within  this  data  should  offer  encouragement  for  experimental  work  at  higher 
impact  velocities  and  into  vastly  different  target  materials. 
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DMENSIONAL  MODEL 


In  an  effort  to  develop  a  logical  set  of  Independent  variables  In 
dimensionally  correct  groupings  agalnat  which  to  test  experimental  data, 
consideration  has  been  given  to  the  trans formation  processes  involved  in 
the  dissipation  of  the  projectile  kinetic  energy  during  a  hypervelocity 
ballistic  impact.  In  the  Third  Quarterly  Report^  the  following  mechanlsaw 
by  which  the  impact  kinetic  energy  may  ultimately  be  dissipated  ware 
discussed.* 

1.  Kinetic  energy  of  the  backsplashed  material 

2.  Energy  required  to  melt,  vaporise,  or  otherwise 
dissociate  projectile  and  target  material 

3.  Energy  emitted  as  electromagnetic  radiation  at  all 
wavelengths 

4.  Energy  transferred  to  atomic,  molecular,  or  granular 
change  such  as  the  energy  of  recrystallisatlon. 

5.  Strain  energy  absorbed  by  the  solid  target  material 

Conservation  of  energy  would  then  require  that  the  impact  kinetic  energy 
be  equal  to  the  sum  of  these  energy  dissipation  terms.  However,  It  became 
apparent  that  possibly  most  of  the  Impact  kinetic  energy  can  be  accounted 
for  in  terms  of  strain  energy  absorbing  parameters. 

Other  than  frictional  effects  which  occur  during  the  process  of  primary 
crater  formation,  the  net  work  done  on  the  target  may  be  determined  largely 
by  the  non- recoverable  or  permanent  deformational  energy  absorbed  by  the 
target  material.  That  is  to  say,  the  energy  dissipated  as  heat,  sound,  elec¬ 
tromagnetic  radiation,  recrystallisatlon,  etc.,  is  essentially  a  result  of 
non- reversible  material  deformation.  Thus,  the  correlation  of  Impact  kinetic 
energy  with  either  mechanical  deformational  work  or  the  combined  chemical, 
mechanical,  acoustical,  thermal,  electrical,  etc.,  energies  of  dissipation 
provides  dual  methods  for  analysing  the  Impact  problem.  Even  the  heat 
generated  by  viscous  effects  on  Impact  cannot  be  dissipated,  except  through 
electromagnetic  radiation,  during  the  very  short  period  of  time  over  which 
the  primary  crater  is  formed.  Since  the  amount  of  electromagnetic  radia¬ 
tion  which  can  take  place  during  crater  formation  is  small  *,  this  heat 

*  A  brief  discussion  of  the  first  four  terms  appears  In  Appendix  A. 
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merely  contribute!  to  the  environment  under  which  the  terget  meterlel 
le  deformed.  The  mejorlty  of  the  beckspleehed  meterlel  roey  be  due  to 
elestlc  recovery  which  occurs  efter  the  projectile  cesses  to  exist  as 
e  ceusltlve  force,  end  la  thus  a  manner  In  which  recoverable  deforms- 
tional  energy  la  dissipated  by  the  target.  Of  course,  there  are  bound 
to  be  other  considerations,  such  es  movement  of  the  entire  target, 
which  are  assumed  to  be  secondary  In  this  treatment. 

The  Idee  that  non- recoverable  strain  energy  accounts  for  the  Impact 
kinetic  energy  before  Its  dissipation  In  other  forms  Is  not  original. 
Goldsmlth^^  points  out  that  two  approaches  have  been  devised  to  account 
for  such  non- reversible  phenomena:  (a)  the  hydrodynamic  theory  and  (b)  the 
theory  of  plastic  flow.  It  la  merely  suggested  here  that  the  parameters 
which  provide  e  measure  of  a  target's  ability  to  "stop"  a  projectile  (and 
In  so  doing,  to  describe  the  cratering  process)  are  appropriate  compression 
and  shear  energy  absorbing  parameters.  This  approach  can  only  be  justified 
by:  (1)  the  fact  that  no  attempt  is  being  made  to  predict  absolute  crater¬ 
ing  (only  observed  variations  in  cratering),  (2)  the  fact  that  static  or 
quasi-static  material  parameters  are  generally  Indicative  of  behavior  under 
vaatly  different  conditions  as  discussed  previously  under  General  Aspects 
of  Target  Behavior.  There  are,  of  course,  discrepancies  Involved  In  the 
use  of  static  or  quasl-statlc  material  parameters,  but  these  difficulties 
are  unavoidable  regerdless  of  the  method  of  enelysls. 

The  problem  of  bridging  the  gap  between  lower  velocity  and  hyper-  .. 
velocity  ballistic  Impact  theory  requires  an  explanation  of  how  the 
projectile  moving  Into  a  target  with  an  Initial  velocity  greater  than 
that  of  any  stress  wave  that  can  be  detached  affects  the  cratering  process. 
One  point  that  seems  evident,  however.  Is  that  In  both  cases,  the  crater¬ 
ing  process  depends  largely  on  the  target's  ability  to  absorb  energy  which 
must  be  related  in  some  way  to  pertinent  impact  and  material  parameters. 

If  the  relations  governing  the  pressures  developed  In  subsonic  and  hyper- 
velocity  Impacts  constitute  the  primary  difference  In  the  two  cases,  it 
Is  shown  In  the  discussion  of  Eqs.  (14,  22,  end  24)  that  the  subsonic  pres¬ 
sures,  determined  from  the  Bernoulli  relation  (Eq.  (14))  and  the  hyper- 
velocity  shock  pressures,  determined  from  one-dlmenslonal  shock  theory 
(Eq.  (22)),  yield  similar  penetration  expressions  In  all  but  possibly 
extremely  high  velocity  Impacts. 

In  order  to  study  the  behavior  of  the  target  during  the  Impact 
process.  It  Is  necessary  to  consider  the  nature  of  the  loading  and  the 
state  of  the  material.  The  ordinary  use  of  the  terms  solid  and  fluid 
seem  to  be  Insufficient  In  describing  the  target  material  under  such  high 
strain  rates  and  pressures,  etc.  The  hydrodynamic  theory  assumes  that 
the  target  material  has  been  stressed  far  beyond  Its  strength  and  there¬ 
fore  behaves  as  a  liquid  with  virtually  no  shear  strength.  However,  since 
extremely  high  pressures  tend  to  congeal  liquids,  the  high  pressures 
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produced  during  Impact  may  be  a  result  of  the  quesl-solld  tet'get  materlel'a 
ability  to  develop  such  high  pressures  following  Impact,  Also,  high  pres¬ 
sures  are  known  to  substantially  Increase  shear  strengths  so  that,  even 
though  some  type  of  materiel  flow  undoubtedly  takes  place,  the  shear  strain 
energy  absorbed  may  bk  a  significant  factor.  Regardless  of  the  magnitude 
of  the  Impact  velocity,  there  must  be  a  certain  period  at  the  end  of  the 
crater  formation  period  In  which  material  strength  plays  a  dominant  role. 

The  compression  and  shear  strain  energy  absorbed  In  the  Impact  process 
depend  generally  on  the  following  factors: 

1.  The  magnitude  and  distribution  of  the  pressures  that 
the  target  material  develops  under  the  Impact  conditions  along 
with  the  rate  at  which  the  high  Initial  pressures  are  attenuated. 

2.  The  compressibility  of  the  target  material  or  Its 
reciprocal,  the  ability  of  the  material  to  resist  compressive 
deformation. 

3.  The  shear  toughness  of  the  target  material  under  the 
Impact  conditions. 

DEVELOIMENT  OF  MODEL 


The  general  expressions  for  compression  and  shear  strain  energy  ab¬ 
sorbed  for  the  case  of  a  simplified  hemispherical  surface  are; 

Compression  Strain  Energy 
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Although  the  extent  of  the  Inelastic  deformation  In  the  target  Is  uncer¬ 
tain,  It  would  appear  to  be  generally  a  function  of  the  target's  ability 
to  absorb  compression  energy.  Thus,  the  upper  limit  (  Tjr  )  for  the  Inte¬ 
grals  might  be  considered  as  some  function  of  the  crater  depth  (constant 
times  crater  depth  (Pc)  for  linear  function). 

A  measure  of  the  unit  volume  change  of  the  target  material,  which 
appears  In  Eq.  (1),  might  be  taken  from  hydrostatic  compressibility  data 
which  has  been  obtained  for  a  few  materials.  An  expression  of  the  form 

A.p  “ 


where  A  and  B  are  constants,  seems  to  yield  an  empirical  fit  with 
existing  data.  It  can  be  noted  In  Table  1  that  the  factor  A  In  Eq.  (3) 
and  the  compressibility  (reciprocal  bulk  modulus)  are  almost  Identical, 
hence 


hf*- 


(4) 

It  la  quite  obvious  that  Eq.  (4)  Is  not  really  applicable  to  very 
high  velocity  Impacts.  Although  the  second  order  term  Is  a  small  cor¬ 
rection  to  the  first  order  term  up  to  a  few  hundred  atmospheres,  at  a 
million  atmospheres  the  two  terms  are  roughly  equal.  Beyond  a  million 
atmospheres,  Eq.  (4)  would  predict  an  Increase  in  volume  with  pressure. 
Dropping  the  second  order  term  to  form 


r 


(5) 


certainly  does  not  fit  the  physical  boundary  conditions  even  at  a  few 
hundred  atmospheres.  Obviously,  the  hugonlot  relationship  developed  for 
a  given  equation  of  state  should  be  used  here,  but  this  Introduces  addi¬ 
tional  unknowns.  It  should  be  re-emphaslsed  that  the  present  model  Is 
only  dimensional  and  that  no  attempt  Is  being  made  for  a  self-consistent 
theory. 


The  manner  In  which  the  pressure  (p)  varies  as  a  function  of  crater 
depth  (r)  Is  not  known.  However,  If  the  assumption  Is  made  that  the  pres¬ 
sure  varies  linearly  from  some  original  pressure  po  at  the  Instant  of 
Impact  (r  ■  0)  to  a  final  pressure  of  zero  when  the  crater  Is  fully 
formed 


(6) 


Then  Eq.  (1)  may  b<s  Integrated  to  yield 

"C,  *■  1^*  (7) 

when  Eqns.  (4)  and  (6)  are  used,  or 
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when  Eqns,  (5)  and  (6)  are  used  ••••  C|  and  C2  are  constants. 

Considering  the  shear  toughness  of  the  crater  material  to  be  the 
product  of  S  and  g  and  the  plastic  shear  strain  In  the  vicinity  of  the 
final  crater  boundary  to  be  a  function  of  this  toughness  factor;  also 
assuming  the  extent  of  the  inelastic  shear  strain  into  the  target  to  be 
soise  linear  function  of  the  crater  depth,  a  solution c£  Eq.  (2)  Is 

(E,V  -CaCSg)^^*  (9) 

where  C3  Is  a  constant. 


Equating  the  Impact  kinetic  energy  to  the  sum  of  Eqns.  (8)  and  (9) 
yields 


or 


(10) 


(11) 


*  Cs(S^)t 

where  C4  and  C3  are  new  constants  and  (hemispherical  crater) 


As  soon  as  the  impacting  materials  begin  to  deform,  the  maximum 
developable  pressure  might  be  considered  as  the  dynamic  fluid  pressure  (or 
stagnation  pressure)  for  the  smaller  density  material. 

(12) 

where  Cg  Is  a  constant  and  Is  the  smaller  of  the  projectile  or  target 
material  density,  ' 


However,  a  better  estimate  of  the  pressure  developed  might  be  obtained 
for  the  case  In  which  both  projectile  and  target  are  regarded  as  Incom¬ 
pressible  fluids.  In  this  case,  from  the  Bernoulli  relation, 

f) -  V*  *  (13) 

where  v  ■  Impact  velocity 

VI  ■  velocity  of  projectile-target  Interface 


From  Eq.  (13), 


or 


Hence, 


V.  .  -jSl - 


(14) 

where  C7  Is  e  constsnt  end  refers  to  a  composite  projectile-target 

density  term  defined  ss 

Interesting  comparisons  are  made  on  Lines  5  and  6  of  Table  1  which  show  the 
variation  of  the  composite  density  term  with  the  projectile  and  target 
densities.  Note  that  when  the  projectile  and  target  densities  are  equal, 
the  composite  density  term  Is  0.25  times  either  one. 

Alternately,  the  Initial  contact  pressure  might  be  determined  by  one- 
dlmenslonal  shock  theory.  The  conservation  of  mass  and  momentum  equations 
are  written  for  a  stationary  shock  front,  where  the  undisturbed  condition 
(ahead  of  the  shock)  and  the  compressed  condition  (behind  the  shock)  are 
denoted  by  the  eubscrlpts  0  and  1,  respectively. 


Conservation  of  Mass 


ft  (u.-Vi)=/^  (4-k) 


C<....r,.tlon  of  ftoeofo.  (u.-V^iv,-t^ 

where  u  Is  the  shock  velocity  and  v  Is  the  particle  velocity. 

Solving  Eq.  (16)  yields  _ 

and 

»V-.  ±J  „e) 

If  the  change  In  compressibility  with  pressure  Is  neglected,  the  shock 
velocity  (u)  can  be  taken  as  the  dllatatlonal  wave  velocity,  or 

U  -  C  (19) 

The  appropriate  expression  for  the  dllatatlonal  wave  velocity  c  for  the 
case  where  the  compressibility  behavior  is  assumed  to  be  linear  and  ^ere 
the  contribution  of  shear  stiffness  G  Is  neglected  (Inelastic  range)  Is 


Solving  Eqs.  (17),  (18)  and  (19),  following  the  assumption  that 
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and  ^9*0  (corresponding  to  the  Initial  iopact  condition)  ylalds, 


for  the  target  *  V,  Ct  /’t 

and  for  the  projectile  a  (  V-  V.;  cys/v 


(20) 

(21) 


where  unsubscrlpted  v  Is  Impact  velocity.  Eliminating  V|  between  Eqs. 
(20)  and  (21)  yields  an  expression  for  the  initial  Interface  pressure 


(22) 


A  comparison  of  the  Bernoulli  pressure  (Eq.  (14))  with  the  shock  pressure 
(Eq.  (22))  In  the  Discussion  Section  shows  that  similar  penetration  ex¬ 
pressions  are  obtained  In  all  but  possibly  extremely  high  velocity  laipacts. 


Substituting  Eq.  (14)  for  In  Eq.  (11)  yields 

.  / _ 


Also  substituting  Eq.  (14)  for  ^  In  Eq.  (10)  yields 

p.  .  _ ^ _ 

The  use  of  Eq.  (7)  rather  than  Eq.  (8)  in  the  development  would  result 

J  ^ 


(23) 

(24) 


In 


v< 


and 


(23*) 

A. _ _ 

The  use  of  Eq.  (22)  rather  than  Eq.  (14)  In  Eqns.  (11)  and  (10)  gives 

*  (25) 


- 


and 


Tt  ‘^10  (sjX 

ni?  MonvT  mMC*rAiinK  ' 


(26) 


AMALYSIS  OF  MODEL  CONSTANTS 


A  preliminary  anlysis  of  the  dimensional  model  constants  was  made  using 
the  Utahr  data.  These  data  considered  a  variety  of  target  materials  although 
only  over  a  limited  range  of  impact  velocities  around  2  Km/sec.  The  slope 
of  the  experimental  crater  voluise /kinetic  energy  curve  (^A^)  from  reference 
11  was  used  as  an  average  value  for  each  material  In  Table'  1.  The  model  In 
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Table  1.  Miscellaneous  Calculations 


Steel  Projectile*  -  Target  Material  - 4140  Steel  Copper  Zinc 


V 

1,  From  "best  fit"  curves  (Utah,  p.l4)-=£  - 

0.20  X  10-1“ 

0.60 

0.65 

2.  (Line  1.)  x  (0.88  x  10^®  dyne- cm)  =  "Vc*  ~ 

0.176  crfi^ 

0.53 

0.57 

3.  = 

7.90  gm/cpi'^ 

8.95 

7.14 

«.  (See  E,.  15  ),  = 

1.96  <^yy>/tyr)^ 

2.11 

1.86 

5.  =.  (Line  4.)//0^  r 

0.25  dimensionless 

0.24 

0.26 

0.25  dimensionless 

0.27 

0.24 

^  (The  steel  value  was  given  for  iron,  >/  _ 
but  value  also  good  for  steel) 

1.67  X  10^^  dijne/cm^ 

1.33 

0.57 

8.  Compressibility;  ^/(Line  7.)=  “ 

0.600  X  10  cni^dijne 

0.752 

1.76 

9.  (Line  8.)  x  (980)  -  V ~ 

0.589  X  10‘^  C.r>iy3>>7 

0.737 

1.722 

]0.  Parameters  for  Eq.  (3  ), 

0.584  X  lO'^  ttry<y*\ 

0.732 

u.  J  = 

0.040  X  10'^'^  Qtri^/<yr\^ 

0.27 

12,  St  = 

55.0  X  10®  Jt^ne/cD?*" 

15.8 

13,1 

13.  Values  for  g  could  not  be  Qa  - 

found.  Tensile  elongation  values  ^ 

0.17  dimensionless 

0.20 

0.15 

were  used.  Both  provide  a  relative 
measure  of  ductility. 

Calculations ; 

14  X  1 

0.074  dimensionless 

0.322 

0.63: 

St 

0.1100  dimensionless 

0.0948 

0.08: 

"■  ^  W 

17,  (Bee  Eq.E7),  Const. =  y  ^  = 

0.0187  dimensionless 

0.0190 

o.oi: 

9.1  dimensionless 

10.5 

11.7 

18.  (See  Eq,^8),  Const.  =  “ 

5.4  dimensionless 

5.3 

7.8 

*  for  /^*  7*80  3  =  0.056V  c»>i^ 

•  Maximum  Deviation  From  Mean  Value. 


Table  1.  Miscellaneous  Calculations 


Copper 

Magn, 

Alum. 

Lead 

Reference 

iterial  ■ 

V 

P.14)l  . 

0-20  *  10 

0.60 

0,65 

0.84 

1.67 

4.25 

Utah 

■-cm)=  = 

0.176  ctn’ 

0.53 

0,57 

0.74 

1.47 

3.74 

Calculated 

A  = 

7.90  gi^y/cm"* 

8.95 

7,14 

1.74 

2.70 

11.4 

Physics  >e 
Handbook  ^ 

1.96 

2.11 

1,86 

0.776 

1.07 

2.33 

Calculated 

0,25  dimensionless 

0.24 

0.26 

0.44 

0.40 

0.20 

Calculated 

= 

0.25  dimensionless 

0.27 

0.24 

0,10 

0.14 

0.30 

Calculated 

jr  iron,  </  _ 
»1)  ^  ~ 

1.67  X  10^^  di^ne^m^ 

1.33 

0.57 

0. 34 

0.75 

0.44 

Kent  9^ 

Handbook 

•  )=Vk,  = 

0.600  X  10  CW^di^ne 

0.752 

1.76 

2.94 

1.33 

2.28 

Calculated 

=  Vk,  = 

0.589  X  lO'^  CmVgW 

0.737 

1.722 

2.89 

1,305 

2.23 

Calculated 

A,  - 

0,584  X  lO'^  Clwygyn 
0.040  X  10' 1’ 

0.,732 

1.334 

2.37 

Nedai, 

t 

Reinhart  jg 
&  Pearson 

0^ 

ti¬ 

ll 

0.27 

0.35 

1.73 

s,  - 

55.0  X  10®  c/t^ne/cm^ 

15.8 

13.1 

13.0 

6.5 

1.2 

Utah 

lues 

0.17  dimensionless 

0.20 

0.15 

0.16 

0,12 

0.30 

Metals  29 
Handbook 

elative 

(Except  ANC-5 
for  Steel) 

0.074  dimensionless 

0.322 

0.633 

0.233 

0.406 

8.419 

= 

0.1100  dimensionless 

0,0948 

0.0852 

0.1092 

0.1086 

0.0510 

= 

0.0187  dimensionless 

0.0190 

0.0128 

0.0175 

0.0130 

0.0153 

9,1  dimensionless 

10.5 

11.7 

9.1 

9.2 

19.6 

707,® 

(s,)*  ' 

5,4  dimensionless 

5.3 

7,8 

5.7 

7.7 

6.5 

227.® 

Vf,  -  0.056V  j  V  =  2.  ^m/stc  ^  Fp  -  ^  =  0,88  ~  cno. 

in  Value, 
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(23) 


the  form  of  Eq,  (23)  was  used. 

Vc  . _ ! _ 

The  various  values  of  the  quantities  in  Eq,  (23)  are  shown  in  Table  1  along 
with  the  appropriate  reference  comments.  One  shortcoming  (noted  on  Line 
13  of  Table  1)  is  the  necessity  of  using  tensile  elongation  values  instead 
of  ultimate  shear  strains  g.  This  is  not  too  unreasonable,  however,  since 
both  provide  a  relative  measure  of  ductility.  The  ultimate  shear  strains 
for  the  different  target  materials  could  not  be  found. 


The  Utah  report  “  showed  an  approximate  correlation  of  volume/ energy 


with  shear  strength(Sj.)  for  all  target  materials  tested  except  lead. 


Jlc 


/ 


Cohst.  CsX 


(27) 


These  constants  are  repeated  on  Line  17  of  Table  1  where  maximum  devia¬ 
tion  from  the  mean  value  is  707..  Since  lead  is  the  most  ductile  material 
of  the  group,  the  possibility  of  a  closer  correlation  of  volume/energy 
with  shear  toughness  is  suggested  as 

=  cJt. 

These  constants  are  shown  on  Line  18  of  Table  1  where  the  maximum  deviation 
from  the  mean  value  is  only  227.,  using  the  tensile  elongation  parameter 
for  the  ultimate  shear  strain  g^. 


Rewriting  Eq.  (23), 


Vc  I 

ICt  V'-t 


or,  referring  to  Table  1, 


(29) 


(line  14)  a^  +  (line  16)  a2  =  1  (30) 

Approximate  values  of  a^^  =  0.02  and  a2  =  60  were  found  by  iteration  using 
the  Utah  data.  When  substituted  into  Eq.  (30),  these  constants  yield  a 
maximum  deviation  from  the  mean  of  207.  for  the  different  materials,  as 
shown  below. 


steel : 

.074 

copper; 

.322 

zinc: 

.633 

magnesium: 

.233 

aluminum: 

.406 

lead; 

8.419 

aj^  +  0.0187  a2  =  1.1 
a^  +  , 0190  32  =  1.1 

a^  +  .0128  a2  =0.8 

ai  +  .0175  32  =  1.1 

aj,  4-  .0130  32  =  0.8 

aj^  +  .0153  32  =  1.1 
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It  Is  Interesting  to  note  In  the  above  relations  that  the  very  much  higher 
coefficient  of  a^  for  lead  may  Indicate  the  relative  liq>ortance  of  the 
density  term  when  the  Impact  velocity  is  in  the  vicinity  of  the  target 
dllatatlonal  wave  velocity  (the  case  only  for  lead  here).  The  maximum  Is^ 
pact  velocities  were  about  2  km/sec. 

OTHER  MOraiLS 

This  section  Is  adapted  from  the  excellent  paper  by  Herrmann  and  Jones 
for  purposes  of  comparison  and  discussion.  By  referring  to  an  arbitrary 
target  realstance  force  F,  the  work-klnetlc  energy  expraaslon  can  be  written 


V»  -FcLPc 


which  can  be  Integrated  to  yield  a  penetration  expression  for  a  given 
definition  of  F.  This  approach,  however,  does  not  Intuitively  Include  any 
work  done  on  the  target  In  the  form  of  shear  deformation. 

Virtually  any  power  expression  of  the  form 

(32) 

can  be  obtained  from  Eq.  (31)  by  assuming  the  resistance  force  to  be  of  the 
form  /  ^  1.  ^ 

Fs  C(  (Pc  j 

where  ji  ^  yields  a  logarithmic  relation).  For  example.  If 

the  resistance  force  Is  assumed  to  be  proportional  to  penetration,  as 

Cjt  Pi 

or  proportional  to  all  four  parameters,  as 

Eq.  (31)  reduces  In  both  cases  to 

^.C,  (33) 

when  the  constant  of  Integration  Is  assumed  zero. 


If  Is  defined  as  some  target  material  strength  parameter,  one 
possibility  that  provides  an  approximate  data  fit  Is 

(3 

where  the  strength  parameter  H^  is  the  Brlnell  hardness.  Eq.  (33)  then 
becomes  ^  #-  a. 

.fl 

nr  ^  u  Urn  (3 


which  can  be  written  In  dimensionless  form  as 
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Bq.  (36)  can  also  be  obtained  by  equating  the  volune/energy  ratio  to 
•nd  assuming  a  hemispherical  crater. 

A  closer  fit  than  Eq.  (36)  has  been  found  by  adjusting  the  power  of  the 
density  ratio  to  _  ^  8/.  .  7.  (A 


\  Ht  ' 


This  expression  Involving  both  density  and  strength  terms  was  shown  by 
Herrmann  and  Jones  to  provide  the  best  power  function  fit  for  a  variety  of 
target  and  projectile  materials  (a^  about  0.36),  although  still  over  a 
somewhat  limited  range  of  Impact  velocities. 

If  now  the  resistance  force  F  Is  assumed  to  be  dependent  on  the  Inertia 
forces  exerted  by  the  target, 

dF.l<lm)a  -.  (Area 

/°t  iv 

•nd  V  d-V 


kCtVf  “Vt 


which  when  substituted  Into  Eq.  (31)  yields  upon  Integration 

(40) 

If  Cg  Is  Interpreted  as  a  material  strength  parameter,  Independent  of 
velocity,  another  expression  similar  to  Eq.  (39)  might  be  written  ea 

r  =  Cl  Vp +  £?,<,  Ht)  (41) 

which  leads.  Instead  of  Eq.  (40),  to 

f  /?  1/  /4  ]*C  II  (42) 

Assuming  the  constant  ^11  to  be  zero, 

(^)  ^ 

Eq.  (43)  was  written  by  Herrmann  and  Jones  as  - 

=Cil  )]  <«) 

where  multiplier  (a  constant  for  a  given  material)  was  delated.  From 
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Eq.  (44),  a  least  squares  fit  of  existing  experlmantal  data  to  an  equation 

of  the  form  o  i  .  1 

VC  _  ^  s  r .  /  \l 


=C  K  1a  [1  + J  )]  („) 

a  more  appropriate  expression  than  Eq.  (44)  to  be 

le  approximated  by  the  expression 

:lty  range  limited  to  medium  and  moderately  high  Impact  velo- 


for  a  velocity  range  limited  to  medium  and  moderately  high  Impact  velo¬ 
cities.  Eq.  (47)  compares  very  closely  with  Eq.  (37). 

The  logarithmic  expression  Involving  both  density  and  strength  terms 
has  been  shown  by  Herrmann  and  Jones  to  fit  experimental  data  over  the 
entire  experimental  range,  but  only  when  appropriate  values  for  the  con¬ 
stants  are  used  for  given  target  and  projectile  materials.  These  constants 
were  found  to  vary  widely  between  materials. 


DISCUSSION  AND  COMPARISON  WITH  EEPERIMENTAL  DATA 

The  dimensional  model  was  checked  against  the  experimental  data  using 
the  second  "low"  velocity  distribution  (Fig.  2)*  and  the  "high"  velocity 
distribution  (Fig.  3)*.  The  least  squares  values  of  the  constants  (a^) 

In  Eqs.  (24),  (25),  and  (26)  were  evaluated  on  the  digital  computer.  In 
addition,  the  multiple  correlation  coefficients  (mcc)  between  the  Inde¬ 
pendent  and  dependent  variables  taken  as  a  whole  were  computed  for  each 
equation.  The  results  were  as  follows: 


Using  the  second  "low"  velocity  distribution 


for  Eq.  (24) 

a3  -  -1.02  a4  -  24.1 

mcc  ■  -0.057 

for  Eq.  (25) 

a?  -  -1.06  x  10-3  ag  -  1.92 

mcc  ■  -0.067 

for  Eq.  (26) 

Sj  -  -0.115  SjQ-  27.5 

mcc  -  -0.080 

the  "high"  velocity  distribution 

for  Eq.  (24) 

a3  -  2.97  X  10*^  a4  -  10.2 

mcc  ■  0.019 

for  Eq,  (25) 

ay  -  6,07  x  10-4  ag  -  1.92 

mcc  ■  0,339 

for  Eq.  (26) 

Sg  -  2.86  X  10-5  ajo**  10,3 

mcc  ••  0.003 

^Appendix  B 
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The  are  not  dlmenalonless  and  their  units  may  be  determined  from  the 
units  of  the  independent  variables  as  listed  in  Table  XIV*. 


The  multiple  correlation  coefficients  are  all  quite  small  indicating 
that  none  of  the  equations  are  doing  a  good  Job  of  predicting  penetration 
over  either  of  the  velocity  ranges.  The  negative  correlation  coefficients 
merely  indicate  that  penetration  goes  down  as  at  least  one  of  the  inde> 
pendent  variables  (probably  velocity)  goes  up.  However,  since  the  corre¬ 
lation  coefficients  are  so  small,  it  is  neither  surprising  nor  significant 
that  some  are  negative. 


The  least  squares  constants  ai  when  substituted  into  their  proper 
equations,  make  the  first  term  in  the  denominator  negligible  in  comparison 
with  the  second,  especially  for  the  "high"  velocity  data.  If,  for  example, 
the  first  term  in  the  denominator  of  Eq.  (24)  is  neglected,  the  resulting 
equation 


i 


v/>  <^1,  ^  ^^94 


^5 


compares  favorably  with  the  exponents  given  on  the  fifth  line  of  Table  XIII 
in  all  but  the  exponent  of  projectile  density. 


At  the  same  time  an  attempt  was  made  to  obtain  a  least  squares  value 
for  the  constant  cj  in  the  form  of  the  Herrmann  and  Jones^  equation  given 
by  Eq.  (36).  This  attempt  met  with  complete  feilure  since  the  constant 
seemed  to  vary  randomly  as  more  and  more  experimental  data  points  were 
added.  Multiple  correlation  coefficients  calculated  using  a  C3  of  the 
order  of  magnitude  predicted  by  the  computer  were  very  low.  Both  Che 
Dimensional  Model  and  the  Herrmann  arid  Jones  equation  suffer  by  compari¬ 
son  with  the  Empirical  Model  in  that  they  do  not  leave  enough  constants 
to  be  eveluated  by  least  squares. 


*Appendix  B 
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SlttMABY  OF  RESULTS  AND  CONCLUSIONS 


The  following  la  a  brief  auomery  of  the  areaa  Inveatlgated  end  the 
reaulta  obtained  In  the  atudy  of  target  penetration  prediction  by  high 
■peed  and  ultra  high  apeed  balllatlc  Impact. 

1.  Exlatlng  experimental  data  relative  to  balllatlc  Impact 
waa  evaluated  on  a  atatlatlcal  baala  ualng  a  digital  computer. 

2.  The  llmltatlonaof  preaent  data  and  experimental  proce- 
durea  were  dlacuaaed  and  areaa  for  future  experimental  work 
were  delineated. 

3.  On  the  baala  of  the  atatlatlcal  atudy,  engineering  formu¬ 
lae  were  developed  which  will  predict  with  an  accuracy  of  90%, 

(a)  the  depth  of  penetration 

(b)  the  volume  of  the  crater 

In  aeml-lnflnlte  targeta  of  common  metala  over  a  wide  range  of 
Impact  velocltlea.  The  application  of  theae  equations  In  pre¬ 
dicting  the  alze  and  velocity  of  projectiles  Impacting  orbital 
structures  la  also  Indicated. 

4.  "Strength"  parameters  In  general  (at  least  for  ultimate 
strength,  shear  strength,  yield  strength  and  Brlnell  hardness 
which  were  tested  In  the  study)  were  found  to  be  statistically 
Indistinguishable  In  empirically  explaining  penetration  behavior. 

5.  A  general  study  of  material  behavior  under  static  and 
relatively  low  strain- rate  loading  and  under  conditions  of  high 
■peed  ballistic  Impact  was  made.  Qualified  Justification  was 
given  for  the  use  of  static  or  quasl-statlc  material  parameters 
In  empirical  and  seml-emplrlcal  analyses  of  high  speed  Impacts. 

6.  A  dimensional  penetration  model  was  developed  from  a 
work-energy  consideration  which  suggests  that  the  nonrecoverable 
target  compression  and  shear  strain  energies  account  for  most  of 
the  kinetic  energy  of  the  projectile. 

7.  For  limited  ranges  of  impact  velocity,  the  dimensional  model. 
Involving  both  compressibility  and  shear  toughness  parameters,  was 
shown  to  yield  consistent  results.  However,  when  extended  to  In¬ 
clude  the  entire  range  of  "experimental"  velocities,  and  for  a 
variety  of  projectile  and  target  materials,  the  predominant 
parameter  proved  to  be  the  shear-strain  energy  absorbing  (shear 
toughness)  parameter,  for  which  the  correlations  obtained  are 
believed  to  be  significant. 
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APPENDIX  A  -  SOKE  METHODS  OF  tNERGY  DISSIPATION 


The  following  are  some  general  conmente  on  four  way*  by  which  the 
Impact  kinetic  energy  may  eventually  be  dissipated  by  the  target. 

1.  Kinetic  energy  of  the  material  backsplashed  during 
cavity  formation 

Eg  B  (1/2  mv^)  backsplash 

Various  opinions  have  been  expressed  as  to  the 
relative  magnitude  of  the  kinetic  energy  of  this 
ejected  material.  It  does  appear  that  the  mass  of 
the  backsplash  material  Is  usually  much  greater  than 
that  of  the  original  projectile.  Although  the  Utah 
group^^  has  measured  very  high  velocities  (up  to  20 
times  the  Impact  velocity  for  mlcron-slzed  spray  par¬ 
ticles),  the  average  velocity  of  the  backsplashed 
material  must  be  quite  small.  This  Is  verified  by 
photographs  of  Impacts  which  show  large  chunks  of 
material  moving  backward  from  the  point  of  Impact 
with  very  small  velocities. 

2.  Energy  required  to  melt,  vaporize  or  otherwise  disso¬ 
ciate  projectile  and  target  material  during  the  crater 
formation.  It  Is  considered  of  little  consequence  here 
that  the  actual  mechanism  of  material  separation  Is  at 
least  In  part  a  result  of  cavitation  following  the 
disappearance  of  Che  projectile.  Bromberg  has  pointed 
out  that  the  terminal  ballistic  data  of  Sunsners  and 
Charters^  can  be  explained  using  the  heat  of  fusion. 


El  "  /®t  ^t  '^p 


If  any  vaporization  takes  place.  It  could  account 
for  a  very  large  portion  of  Che  Impact  kinetic  energy. 

A  simple  calculation  for  an  aluminum  projectile  strik¬ 
ing  an  aluminum  target  at  an  Impact  velocity  of  9  km. /sec. 
shows  that  It  would  Cake  1/3  of  the  Impact  kinetic  energy 
to  vaporize  the  projectile.  If  the  crater  volume  were  57 
times  the  projectile  volume  (not  an  unreasonable  figure) 
the  energy  necessary  to  vaporize  Che  entire  crater  would 
be  over  10  times  the  impact  kinetic  energy.  The  vapor¬ 
ization  term  has  been  eliminated  from  equation  shown 
above  because  It  appears  that  only  a  small  portion  of  the 
projectile  Is  actually  vaporized.  To  create  the  spray  of 
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mlcron-slze.  backsplashed  particles  being  studied 
by  the  Utah*^  group  would  require  at  least  a  portion 
of  the  vaporization  energy,  but  here  again,  the  total 
mass  of  this  material  Is  small. 

3.  The  energy  emitted  by  electromagnetic  radiation  at  all 
wavelengths  due  to  the  heat  generated  at  the  point  of 
Impact.  If  the  absolute  temperature  (T)  of  the  emit¬ 
ting  surface  were  constant  over  an  area  (A)  of  the  target, 
this  radiant  energy  would  be  given  by  the  Stefan-Boltzmann 
law  as 


Ej^  -  e  0-  T*  A  t 

where  (e)  Is  the  emlsslvlty,  (<r)  Is  the  Stefan-Boltz¬ 
mann  constant,  and  t  Is  the  total  time  during  which 
emission  takes  place. 

The  Insignificance  of  this  term  can  be  Illustrated 
by  taking  a  simplified  example.  Consider  an  aluminum 
projectile,  0.64  centimeter  diameter,  Impacting  at  5 
kilometers  per  second.  A  reasonable  ratio  of  crater 
to  projectile  volume  at  this  velocity  would  be  18.  If 
the  crater  rim  area  Is  taken  as  the  emitting  area  and 
the  crater  Is  assumed  to  be  hemispherical,  then  A  =•  3.14 
square  centimeters.  Taking  the  time  over  which  radiation 
takes  place  to  be  the  upper  limit  of  the  time  for  full 
crater  formation  (t  100  microseconds),  the  temperature 
necessary  to  yield  a  ratio  of  Eo/Ep  *  1/1000  would  be 
about  10,000OK. 

4.  Energy  transferred  to  atomic,  molecular,  or  granular 
change  such  as  the  energy  of  recrystallleatlon 

Ea  -  ? 

The  Utah^^  group  has  shown  that  from  12  to  15%  of  the 
energy  required  to  push  a  lead  ball  Into  a  lead  target 
quasl-statlcally  goes  Into  recrystallliatlon.  It  has 
not  yet  been  determined  what  Independent  variables  might 
enter  Into  this  term. 
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Figure  1.  First  "Low”  Velocity  Distribution  of  Target 
Materials  as  a  Function  of  Impact  Velocity 
985  Shots 
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Figure  2  Second  "Low  Velocity"  Distribution  of  Target 
Materials  as  a  Function  of  Velocity  -  742 
Shots 
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Table  VII.  Second  "Low"  Velocity  Data  Evaluation  of  the  Coefficients  in  the  Equation 
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Table  DC,  Second  "Low"  Velocity  Data  Evaluation  of  the  Coefficients  In  the  Equation 
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Figure  3.  "High  Velocity"  Distribution  of  Target  Haterlala 
as  a  Function  of  Velocity  -  266  Shots 
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Table  X.  "High"  Velocity  Data  Evaluation  of  the  Coefficients  in  the  Equation 
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Table  XI.  "High"  Velocity  Data  Evaluation  of  the  Coefficients  kj  in  the  Equation 
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^ble  XII.  "High"  Velocity  Data  Evaluation  of  the  Coefficients  k.  in  the  Equation 
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Table  XIV.  Material  Parameters 
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